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I.  INTRODUCTION 


This  addendum  augments  NOSC  TR  208  (ref.  1 ) and  should  be  read  in  conjunction 
with  that  report.  In  the  report,  a mathematical  model  of  the  wideband  HF  channel  was 
developed  and  used  for  computer  simulation  of  the  channel.  The  purpose  of  this  addendum 
is  to  augment  the  existing  channel  model  by  adding  (1)  a deviative  absorption  term  for  the 
skywave  returns,  and  (2)  an  attenuation  term  for  the  groundwave.  With  these  additions  the 
mathematical  model  will  more  accurately  portray  the  HF  channel  and  its  utility  will  be  con- 
siderably enhanced. 

In  the  following  sections  the  two  additions  are  described  and  a FORTRAN  listing  of 
the  revised  simulator  is  provided.  Notations  and  definitions  used  are,  wherever  possible, 
identical  to  those  appearing  in  the  report. 

II.  DEVIATIVE  ABSORPTION 

The  earlier  channel  model  incorporated  only  three  sources  of  attenuation  for  the 
skywave  returns.  These  were  nondeviative  absorption,  antenna  gain  and  spreading  (l/r^)  loss. 
Several  less-important  sources  of  attenuation,  such  as  polarization  mismatch,  focusing  and 
deviative  absorption,  were  not  included.  While  the  losses  due  to  these  other  sources  are  rela- 
tively small,  one  of  them  - deviative  absorption  - plays  an  important  role  in  signal  analysis 
applications.  Both  the  deviative  absorption  loss  and  the  signal  distortion  are  greater  for  high- 
angle  returns  than  they  are  for  low-angle  returns.  If  deviative  absorption  is  not  included  in 
the  model,  the  high-angle  returns  will  have  an  exaggerated  effect  on  the  overall  received  sig- 
nal, and  this  may  lead  to  a more  pessimistic  evaluation  of  the  channel  than  is  warranted. 

In  what  follows,  a uactable  approximation  for  the  deviative  absorption  loss  is 
described  and  sample  calculations  using  this  approximation  are  presented.  We  begin  with 
some  notation  and  assumptions;  all  of  which  have  been  discussed  in  greater  detail  in  (ref.  1 ). 

The  electron  density  profile  is  assumed  to  have  a parabolic  shape  in  the  E and  F 
regions;  no  provision  has  been  made  for  daytime  splitting  of  the  F region  into  distinct  F j and 
F2  regions.  The  nighttime  electron  density  profile  is  given  by 

f 0 , h<h£-2HE 


hE-2HE<h<hE  + 2HE 


NCh)  = i 


0 


hE  + 2HE<h<hF-2HF 


(1) 


hp  - 2Hp  ^ h < hp  + 2Hp 


^ 0 , hp  + 2Hp  < h 

Here,  NE  is  the  peak  electron  density  (free  electrons  per  m ) of  the  E region  while  hE  is  the 
height  of  the  peak  electron  density  in  the  E region  in  km,  and  2HE  is  the  semithickness  of 
the  E region  parabola  in  km.  Similar  definitions  apply  to  the  F region  parameters  Np,  hp 
and  2Hp. 

1 . NOSC  TR  208,  HF  Channel  Simulator  for  Wideband  Signals  (U),  by  R Lugannani  and  HG  Booker, 
Unclassified,  31  March  1978. 


For  the  daytime  electron  density  profile,  allowance  is  made  for  ionization  between 
the  E and  F regions.  Between  these  regions  the  electron  density  is  assumed  to  be  constant 
and  equal  to  Ne.  For  the  daytime  profile  we  have, 


The  constants  appearing  in  equation  (2)  are  defined  in  the  same  manner  as  they  were  for  the 
nighttime  case,  except  that  their  values  will  be  different  because  of  the  different  heights  and 
ionization  densities  in  the  daytime  ionosphere. 

In  the  absence  of  a magnetic  field,  the  refractive  index  is  given  by. 


2 ’.m  - i r N(>1> 

= /i-(h)  = 1 - K — T- 
f- 


(3) 


where  f is  the  frequency  in  Hz,  N(h)  is  given  by  either  equation  ( 1 ) or  (2),  and. 


K = * s 80.5.  (4) 

4ir  - eQm 

The  penetration  frequency  for  each  region  is  defined  as  the  plasma  frequency  at  the  peak 
electron  density.  For  the  E and  F regions,  these  are  given  by. 


fpE“ANE 
*pF  = ANF 

with  fp£  and  fpp  in  Hz. 

For  vertical  incidence  in  the 
be  used  for  the  deviative  absorption 


(5a) 

(5b) 


absence  of  a magnetic  field,  the  following  expression  will 
loss  (ref.  2) 


, (vert) 
Ldev 


(0  = 


(6) 


2.  Davies.  K,  Ionospheric  Radio  Propagation,  in  National  Bureau  of  Standards  Monograph  80,  U.  S.  Govern- 
ment Printing  Office,  1965. 
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Here,  L^eyrt)(f)  is  the  attenuation  in  decibels  per  hop,  v = ^(h)  is  the  collision  frequency  in 
Hz,  c is  the  velocity  of  light  in  km/sec  and  hr  is  the  height  at  which  p = 0.  In  the  case  of 
oblique  incidence  the  deviative  loss  is  obtained  from  (6)  by  means  of  the  relationship 


^dev  ~ ^dev*^  ” 


1 


L*vert)(f  • cos0) 
cos  d dev 


(7) 


where  0 is  the  ray  angle  measured  from  the  vertical. 

The  integration  in  equation  (6)  is  difficult  to  perform  because  of  the  presence  of  the 
collision  frequency  in  the  integrand.  To  overcome  this  difficulty  the  integrand  will  be  sim- 
plified by  taking  the  collision  frequency  to  be  constant  in  the  E and  F regions  and  allowing 
it  to  decrease  exponentially  between  the  two  regions.  Thus,  for  the  nighttime  collision  fre- 
quency, we  will  use. 


Hh)  = < 


h<hE-2Hp 

hE-:HE<h<hE  + 2HE 

hE  + 2He  < h < hF  - 2Hf 
hp  - 2Hp  < h ^ hp  + 2Hp 
hp  + 2Hp  < h, 

and  for  the  daytime  collision  frequency, 
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,‘/j 


< h < hp  + 2Ht 


hp  + 2Hp  < h 


where 


*/E  = Mhp:) 
vp  = vthp) 


(10a) 

(10b) 


Substituting  equations  (8)  and  (9)  into  (6),  performing  the  integration  and  using  (7) 
we  obtain  for  the  deviative  absorption  loss: 


In  the  foregoing  analysis,  the  effect  of  the  earth’s  magnetic  field  has  been  neglected, 
but  we  take  this  to  be  an  adequate  approximation  for  the  ordinary  wave.  For  the  extraordi- 
nary wave,  the  penetration  frequencies  are  shifted  relative  to  those  for  the  ordinary  wave 
and  it  is  essential  to  take  this  into  account.  As  in  (ref.  1)  the  extraordinary  wave  penetration, 
frequencies  are  obtained  using  the  relationships, 

l'/2 

(16a) 


f(x)  = if 
pE  -> 


:'het  ‘HE 


f(x)  _I 


lpF  2<hf+  I'HF 

Here,  the  superscripts  “0”  and  “x”  denote  the  ordinary  and  extraordinary  waves  respectively, 
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with  f(0J  Riven  by  (5a)  and  f(9?  given  by  (5b).  The  frequencies  fHE  and  fHF  are  the  gyro- 
pE  pr 


frequencies  associated  with  the  earth’s  magnetic  field  at  the  levels  of  maximum  ionization 
density  in  the  E and  F regions  respectively.  For  the  gyrofrequencies  we  will  use  the  approxi- 
mate expressions  (ref.  2), 


fHE  " 0 87  I 


‘HF 


6370  V 

5 Vi 

[1+3  sin-  <t>] 

(17a) 

6370  + hp  / 

6370  V 

3 ->  Vi 

[1+3  sin-  4>] 

(17b) 

6370 + hp/ 

where  fup  and  f^p  ar«  in  MHz,  and  4>  is  the  latitude  expressed  in  magnetic  coordinates.  If 
<p  and  X are  the  respective  geographic  latitude  and  longitude,  and  <f>0  and  XQ  represent  the 
latitude  and  longitude  of  the  north  magnetic  pole  (0O  ^ 78.3°N,  XG  2 69°W),  then  (ref.  2), 


sin  <l>  = sin  0 sin  0 + cos  <t>  * cos  <t>0  ’ cos  (X  - XQ). 


(18) 


We  assume  that,  to  calculate  deviative  absorption  for  the  extraordinary  wave,  we  may 

( X ) ( X ) 

use  equations  (11)  and  ( 1 2),  but  with  F^  and  lpp  substituted  for  the  ordinary  wave  penetra- 
tion frequencies.  Otherwise,  the  calculations  for  deviative  loss  take  the  same  form  as  those 
for  the  ordinary  wave. 

To  evaluate  Ldev,  the  collision  frequency  is  needed  and  for  this  the  following  approxi- 
mation will  be  used. 


y(h)  = expl  1 2.43  - 0.1 773  • (h  - 90)) 

+ K • (3.93  X 105)  • T-3/“  • N(h) 


(19) 


Here,  T = T(h)  is  the  temperature  at  height  h in  degrees  Kelvin  and  K is  given  by  equation  (4). 
The  First  term  on  the  right  side  of  (19)  accounts  for  collisions  between  electrons  and  neutral 
particles;  it  is  based  on  table  9.4  of  (ref.  3).  This  term  is  important  in  the  E region  but  almost 
negligible  in  the  F region.  The  second  term  in  (19)  accounts  for  collisions  between  electrons 
and  ions  and  has  been  taken  from  equation  4.2.10  of  reference  4.  This  term  is  relatively 


3.  Banks.  PM,  and  Kockarts.G,  Aeronomy,  Part  A,  Academic  Press,  1973 

4.  Ratcliffe,  JA.The  Magneto-Ionic  Theory  and  its  Applications  to  the  Ionosphere,  Cambridge  University 
Press.  1959 
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small  in  the  E region,  but  it  is  the  dominant  term  in  the  F region.  For  the  heights  of  interest 
here,  equation  (19)  provides  an  adequate  approximation  to  the  collision  frequency. 

For  the  temperature,  T(h),  we  will  use  the  profile  appearing  in  the  1959  ARDC 
Model  Atmosphere  (ref.  5).  Throughout  the  range  of  interest  in  h,  this  temperature  profile 
will  be  approximated  by  the  polynomial, 

T(h)=  b0  + bjh + b:h2 + b3h3.  (20) 

The  coefficients  have  been  determined  using  a combination  of  least  squares  and  minimax 
fits.  The  results  are  listed  in  table  1. 

The  deviative  absorption  loss  has  been  incorporated  into  the  Fortran  version  of  the 
simulator  and  a revised  listing  of  the  program  is  presented  in  appendix  A.  The  changes  appear 
mainly  in  the  existing  subroutine  “ATTEN”  and  in  the  new  subroutines  “DEVIAT”  and 
“TEMP.” 

A synthetic  ionogram  has  been  generated,  using  the  revised  computer  program.  The 
resulting  curves  are  presented  in  figure  1 and  they  should  be  compared  with  figure  5 of  this 
report  (ref.  1 ).  The  delays  are  identical  for  the  two  figures,  but  the  increased  attenuation  of 
the  high-angle  return  is  readily  evident. 


Table  I.  Coefficients  for  polynomial  approximation 
of  temperature  profile 

Coefficient  100  km  < h < 210  km  210  km  < h < 600  km 


b0 

1 .414  X 10J 

1.424  X 10-’ 

bl 

-4.428  X 101 

3.990  X 10"  2 

b: 

4.217  X 10"' 

-1 .040  X 10"3 

b3 

-1.004  X |0~3 

3.000  X 10'6 

5.  Min/er.  RA,  Champion,  KSW,  and  Pond.  HL,  The  ARDC  Model  Atmosphere:  Air  Force  Surveys  in 
Geophysics,  No.  1 1 5,  1959. 

1 2 34  5 6 7 8 9 10  11 


FREQUENCY  (MHz) 


Figure  I . Computed  ionosphere  returns  for  one.  two,  and  three  hops. 
(Numbers  on  curve  are  attenuation  in  d B.) 
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111.  GROUNDWAVE  ATTENUATION 


This  section  is  concerned  with  the  attenuation  of  a groundwave  that  has  propagated 
across  a rough,  curved,  lossy  ocean  surface.  Our  approach  is  based  on  the  work  of  Barrick 
(ref.  6)  and  represents  an  attempt  to  reduce  his  expressions  to  a tractable  form  for  computa- 
tion. By  providing  for  internal  calculation  of  the  groundwave  attenuation,  the  utility  of  the 
simulator  is  increased;  since,  it  is  no  longer  necessary  to  enter  this  attenuation  as  a separate 
input  whenever  range,  frequency,  or  sea  state  change. 

Following  Barrick.  the  groundwave  loss,  Lg,  will  be  written  as  the  sum  of  two  terms. 
The  first  term,  L*,  \ is  the  basic  transmission  loss  due  to  propagation  across  a smooth, 
curved,  lossy  (4  mhos/meter  conductivity)  surface.  The  second  term,  l1-\  is  the  additional 
transmission  loss  caused  by  surface  roughness.  We  have 


Cl) 


with  all  losses  in  decibels.  Both  the  transmitter  and  receiver  are  assumed  to  be  located  at  the 
ocean  surface.  The  maximum  error  incurred  as  a consequence  of  this  assumption  is  approxi- 
mately 37c  (ref.  6)  and  is  worth  tolerating  in  view  of  the  attendant  simplicity  (viz,  it  is  not 
necessary  to  specify  the  heights  of  the  transmitter  and  receiver). 

For  the  basic  transmission  loss  the  following  expression  will  be  used. 


20  log , o 


+ C,(f)  • R 


C\lf) 


(22) 


where  R is  the  distance  between  transmitter  and  receiver  in  km,  f is  the  frequency  in  MHz, 
and 


(',(0  = 0.78  X ICT4)  • I'2  58  • exp  1-0.40  X 10"8)  • f5]  (23) 

C -.( f ) = 1.83  • r°- 243  (24) 

The  first  term  on  the  right  side  of  (22)  is  the  loss  the  signal  experiences  by  traveling  a dis- 
tance R across  a smooth,  tlat.  perfectly  conducting  surface.  The  second  term  is  a correction 
term  for  the  earth's  curvature  and  for  a lossy  surface;  it  is  based  on  figure  8 of  reference  6. 
The  coefficients  C | (0  and  C'-»(f)  have  been  determined  using  a combination  of  least  squares 
and  minimax  fits  to  give  a good  approximation  to  Barrick's  results  throughout  the  ranges 
and  frequencies  of  interest.  A plot  of  equation  (22)  is  presented  in  figure  2 here,  together 
with  selected  values  from  Barrick.  As  can  be  seen,  the  agreement  is  very  good. 


6.  Bat  telle  Memorial  Institute  Report.  Contract  DAAH01 -70-C-03 1 2.  Theory  of  Groundwave  Propagation 
Across  a Rough  Sea  at  Dekameter  Wavelengths,  by  DE  Barrick,  1970. 
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Figure  2.  Basic  transmission  loss  for  groundwave. 

Determination  of  the  additional  transmission  loss.  Lg“\  requires  a knowledge  of  the 
ocean  wave  spectrum,  which  is  characterized  by  specifying  the  sea  state.  Two  spectra  have 
been  considered  by  Barrick  (ref.  6):  the  Neumann-Pierson  spectrum,  modified  to  account 
for  wind  direction;  and  the  Phillips  spectrum,  which  is  independent  of  wind  direction.  The 
Phillips  spectrum  was  chosen  because,  as  noted  by  Barrick,  it  appears  to  give  better  agree- 
ment with  actual  observations.  Also,  the  Phillips  spectrum  results  in  losses  that  lie  between 
those  obtained  for  the  upwind-downwind  and  crosswind  directions,  using  the  modified 
Neumann-Pierson  spectrum,  and  it  does  not  require  that  wind  direction  be  specified  as  an 
additional  input. 

For  the  additional  transmission  loss,  the  dependence  on  range,  frequency,  and  sea 
state  is  complicated  and  precludes  a simple,  accurate  analytic  approximation.  However, 
interpolation  can  be  used  to  provide  good  accuracy  as  well  as  ease  of  computation  and  it  is 
this  approach  that  will  be  used.  The  interpolation  scheme  is  linear  in  both  range  and  fre- 
quency, with  reference  values  obtained  from  figures  25-31  of  reference  6.  The  reference 
frequencies  used  are  3,  5,  7,  10,  15,  20  and  30  MHz,  while  the  reference  ranges  used  are  25, 
50,  100,  200,  300,  400,  500,  600,  700,  800,  900  and  1000  km.  With  these  reference  points 
the  error  attributable  to  the  interpolation  is  negligible  for  the  ranges  and  frequencies  of 
interest. 

A Fortran  version  of  the  calculations  necessary  to  determine  the  groundwave  loss 
has  been  incorporated  into  the  simulator;  a listing  of  this  program  is  presented  in  appendix  A 
of  this  addendum.  The  changes  consist  of  a new  subroutine  “GNDATN”  and  a statement  in 
the  main  program  that  calls  this  subroutine. 
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IV.  CONCLUSION 

Two  additions  have  been  made  to  the  HF  channel  model  and  the  associated  simulator. 
These  are: 

( 1 ) Inclusion  of  deviative  absorption  loss  for  the  skywave  returns,  and 

(2)  Automatic  computation  of  the  groundwave  attenuation. 

The  first  of  these  improves  the  accuracy  and  utility  of  the  mathematical  model,  while  the 
second  eliminates  the  need  to  enter  a new  value  for  groundwave  attenuation  whenever  range, 
frequency,  or  sea  state  change. 


GLOSSARY  OF  NEW  MATHEMATICAL  TERMS* 

Coefficients  used  in  the  polynomial  approximation  of  the  temperature  profile 
Parameters  used  in  the  approximation  of  L^* ) 

Deviative  absorption  loss,  oblique  incidence 
Deviative  absorption  loss,  vertical  incidence 
Groundwave  loss 

Basic  transmission  loss  component  of  Lg 
Additional  transmission  loss  component  of  Lg 


T = T(h)  Temperature 

v(h)  Collision  frequency 


Collision  frequency  at  height  hg 
Collision  frequency  at  height  hp 


*This  glossary  contains  only  those  symbols  that  are  introduced  for  the  first  time  in  this  addendum.  Other 
symbols  used  in  the  text  are  listed  in  the  glossary  of  the  report  (ref.  1). 
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